Vortex duality: watching the dual side with order propagators. 
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In condensed matter physics, Kramers- Wannier duality implies that the state disordered by quan- 
tum fluctuations or temperature actually corresponds with an ordered state formed from the topo- 
logical excitations of the 'original' ordered state. At first sight it might appear to be impossible 
to observe this dual order using means associated with the original order. Although true for Ising 
models, we consider in this paper the well known vortex duality, in particular in the quantum in- 
terpretation in 2-l-lD where it is associated with the quantum phase transition from a superfluid 
to a Bose Mott insulator. Here, the disordered Mott insulating state is at the same time a dual 
superconductor corresponding with a Bose condensate of vortices. We present a simple formalism 
making it possible to compute the velocity propagator associated with the superfluid in terms of the 
degrees of freedom of the dual theory. The Mott insulator is characterized by a doublet of massive 
modes and we demonstrate that one of these modes is nothing else than the longitudinal photon 
(gauged second sound) of the dual superconductor. For increasing momenta, the system rediscovers 
the original order, and the effect on the velocity correlator is that the longitudinal photon looses its 
pole strength. The quantum critical regime as probed by the velocity correlator is most interesting. 
We demonstrate that at infinite wavelength the continua of critical modes associated with second 
sound and the dual longitudinal photon are indistinguishable. However, at flnite momenta they be- 
have differently, tracking the weight reshuffling found in the quasiparticle spectrum of the disorder 
state closely. 

PACS numbers: 75.10.Jm, 75.40.Gb, 05.30.Jp 



I. INTRODUCTION 



The notion of Kramers- Wannier duality^ has been 
around for a long time in statistical physics and field 
theory but one can still wonder if it is scope is fully ap- 
preciated. Especially in condensed matter physics it ap- 
pears to be much more than a mathematical convenience. 
One can view it instead as a physical 'relativity' princi- 
ple associated with order. Order and disorder have no 
objective meaning but just depend on the viewpoint of 
the observer. For instance, according to the classic 2D 
Ising model duality^, an observer equipped with machin- 
ery measuring two point correlators of the order degrees 
of freedom (the Ising spins) will be convinced that the 
low temperature state shows long range order while the 
high temperature state is just an entropy dominated fea- 
tureless entity. On the other hand, the Kramers- Wannier 
duality demonstrates that the high temperature state in 
fact corresponds with an ordered state, a condensate, 
formed from the topological excitations (domain walls) 
associated with the low temperature order. An exper- 
imentalist probing the system with a machine sensitive 
to the domain wall order would be of the opinion that 
the high temperature state is ordered, while the low tem- 
perature state is entropy dominated. When duality is in 
charge, disorder is order in disguise and one might think 
that this camouflage act is perfect. The dual order is 
carried by the topological excitations of the direct order. 
In the continuum limit it takes an infinity of operations 
involving the order degrees of freedom to probe the topo- 
logical excitations and this is beyond the capacity of any 



machine builder. Henceforth, it appears that it is fun- 
damentally impossible for an observer which can only 
employ order degrees of freedom to directly measure the 
order associated with the dual side with the consequence 
that the 'order-experimentalist' can only perceive dual 
order as disorder. 

Here we will demonstrate that the above is too strong 
a statement. This 'duality censorship' seems absolute for 
the special case of the Ising model in 2D. However, in 
a way its scalar order parameter structure is too sim- 
ple. Here we will focus on a more representative ex- 
ample: the 3D XY model which might be alternatively 
interpreted as the Bose-Hubbard model in 2-l-lD at zero 
chemical potential^, or either as the Abelian-Higgs model 
in 2-l-lD of high energy physics at T = 0. It is charac- 
terized by the well-known Abelian-Higgs'* or vortex du- 
ality which maps the global XY model on [/(1)/C/(1) 
gauge theory. In the Bose-Hubbard interpretation, the 
quantum disordered neutral superfluid corresponds with 
a dual Meissner phase characterized by Bose condensed 
vortex-particles. This incompressible state corresponds 
physically with the Bose Mott-insulator— . 

The excitations of the 'dual side' are the Higgs (am- 
plitude) mode and massive photons of the dual super- 
conductor. Surely, the Higgs mode is subjected to dual 
censorship for the same reasons as found in the Ising 
model, but the photons are a different story. As is well 
known, the Goldstone mode (second sound) of the su- 
perfluid turns into a doublet of massive excitations in 
the Mott-insulator corresponding with the propagating 
unoccupied- and doubly occupied states in the Mott 
state. As we will show here, these actually correspond 



with linear combinations of the degenerate transversal- 
and longitudinal photons of the dual superconductor. 
The latter is of course the usual extra gauge mode as- 
sociated with the presence of the dual phase order, and 
one can say that the dual censorship does not prohibit 
the dual phase order to manifest itself on the order side. 

Our working horse is a simple expression relating the 
second sound propagators to the dual photons which just 
follows from the Legendre transformation (section III). In 
combination with the essentially complete understand- 
ing of the physics resting on the dual order, we are able 
to describe with little effort some features of the order 
propagators in the disordered state which are to the best 
of our knowledge not recognized. On the gaussian level 
(section IV, V) we obtain the outcome sketched in Fig. 
^p: the single Goldstone of the ordered state (Fig. 
turns into a massive doublet in the disordered state. At 
zero momentum the two poles of the velocity propagator 
are indistinguishable but one of them looses its strength 
when momentum is exceeding the inverse London length 
of the dual superconductor. This makes sense because at 
shorter lengths order is re-emerging and the simple Gold- 
stone spectrum of the ordered state should be recovered. 

The Abelian Higgs model in 2-l-lD is below its upper 
critical dimension and its critical state is strongly inter- 
acting. Resting on the complete description due to Hove 
and Sudb0^ of this critical state in the dual language, we 
will derive the critical velocity-velocity propagators (sec- 
tion VI) with the surprising outcome that its transversal- 
and longitudinal components appear to be quite different 
although governed by the same anomalous dimension, 
again reflecting the rather different status of the 'order' 
(transversal) and 'disorder' (longitudinal) photons when 
measured through velocity correlations. Before getting 
deeper into this, let us start with some simple consider- 
ations regarding the mode counting. 

II. A SIMPLE COUNTING ARGUMENT 

The system of interest is the well-known Bose-Hubbard 
model in 2-l-lD at vanishing chemical potential, written 
in phase-number representation as^, 

defined on a 2D bipartite lattice; Ui and (jfi are the 
number- and phase operators on site i, satisfying the 
commutation relation [ni,(j)j] — iSij. The first- and sec- 
ond term in Eq. ^ represents the charging- and Joseph- 
son energy respectively. When the coupling constant 
g = 2/{JC) is small the Josephson energy will dominate 
and the phase is ordered at zero temperature, while the 
excitation spectrum consists of a single Goldstone mode 
(phase mode or second sound) shown in Fig. c. On 
the other hand, when g is large the phase is quantum 
disordered and number condenses such that = mod- 
ulo local fluctuations, signaling the Mott-insulator. Of 
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FIG. 1: The excitations in the weak/strong coupling limits 
of the Bose Hubbard model at zero chemical potential: The 
Goldstone boson (second sound) with linear dispersion (c) as- 
sociated with the superfluid (phase ordered) state (a) at weak 
coupling. In the strong coupling limit (b) a doublet of mas- 
sive 'excitons' are found with gap A (d) corresponding with 
propagating unoccupied and doubly occupied sites, which can 
be alternatively understood at g — + as the ±1 angular mo- 
mentum eigenstates of a 0(2) quantum rotor. 



central interest is the excitation spectrum of the Mott- 
insulator. In the rotor language^, the ground state is the 
angular momentum singlet while the lowest lying excita- 
tions consist of a doublet of propagating M = ±1 modes 
characterized by a zero- momentum mass gap (Fig. ^). 
In the Bose-Hubbard interpretation these have a simple 
interpretation in the strong coupling limit (g — > oo) as 
bosons added (A/ = +1) or removed {M = — 1) from the 
charge-commensurate state (Fig. ^), while their delo- 
calization in the lattice produces an identical dispersion 
due to the charge conjugation symmetry of the model 
Eq. (P). 

This is of course overly well known, but one might be 
at first a bit puzzled about how to relate these excitations 
to the well established notion that the phase disordered 
state does correspond with a dual (vortex) superconduc- 
tor. After all, the doublet of massive modes are fully 
protected elementary excitations and they should arise 
regardless the way one wants to describe the system. 

In the scaling limit, sufficiently close to the quantum 
phase transition where a continuum field-theoretic de- 
scription applies, phase dynamics might as well be de- 
scribed in terms of the dual 'disorder field theory'. This 
turns out to be just the Ginzburg-Landau- Wilson theory 
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of a relativistic 2+lD U{1) superconductor—, 

iTO2|^y|2+w|*v|4 (2) 

The Higgs field ^'y describes the Bose condensate of 
vortices, i.e. the tangle of vortex worldlines. In 2+ ID 
vortex 'particles' are in a precise sense indistinguish- 
able from electromagnetically charged particles. The 
long range interactions between vortices mediated by 
the phase condensate can be described in terms of non- 
compact U{1) gauge fields A^, U{1) appearing in the 
connections for the matter field and in the field strength 
F^i, = d^Au — d^Afj^. Of crucial importance for what 
comes, Eg.© is fully relativistic. The phase-dynamics 
problem is characterized by a single (spin-wave) velocity 
which we take to be one. Although in dual representa- 
tion we are dealing with two separate fields \E'y and , 
both matter- and gauge fields are governed by the same 
velocity. 

The ordered- and disordered phase correspond with the 
Coulomb- and Meissner phase, respectively, of the theory 
Eq. (O . Obviously, the observable consequences of the 
dual theory Eq. Q and the original theory Eq. ^ have 
to be the same and for this to be the case, it is a necessary 
condition that the mode content of both theories are the 
same. 

Let us first consider the superfluid phase/dual 
Coulomb phase. Since phase is condensed, the sys- 
tem should be characterized by a single Goldstone bo- 
son (spin- wave/second sound). How to count this on 
the dual side? The dual Lagrangian is just the 2+ ID 
non-compact U(l) Maxwell Lagrangian This is 

characterized by three vector potentials Ar^A^yAy and 
one gauge constraint. Accordingly, the dual theory is 
characterized by two physical degrees of freedom. This 
is confusing at first, but as we will discuss in some de- 
tail in the next section, it makes a perfect sense. One 
photon is dynamical and transverse ( At = ^^^ ) and 
is just the Goldstone mode in the dual language. The 
other physical photon is the temporal one A-r , describ- 
ing the instantaneous ('Coulomb') interaction between 
static vortex sources: in the Coulomb gauge the longitu- 
dinal photon {Al = ^^) drops out as we subject it to 
the Coulomb gauge fix V • A = 0. Hence, we find the 
correct mode content: the gauge field description is just 
more efficient than the phase description, because the in- 
stantaneous vortex-vortex interactions appear explicitly, 
while they have to be constructed by hand in terms of 
the phase degrees of freedom. 

The full powers of the dual gauge theory unfold in the 
phase-disordered state. In terms of phase, all one has 
is a strong coupling expansion in the Hamiltonian lan- 
guage of Eq. ^ we already alluded to, leading to the 
conclusion that at least in the long wavelength limit one 
is dealing with a twofold degenerate massive 'exciton'. 
What is the dual gauge theory telling us? As will be- 
come clear later, the vortex amplitude 'Higgs boson' \'^v\ 



is subject of dual censorship, and what remains are the 
original photons and the dual phase degree of free- 
dom (pv defined through '^v — |^'y|e'"^^. If the system 
would be non-relativistic such that cv/c —>■ oo, where cy 
and c are the vortex condensate phase velocity and spin 
wave velocity, respectively, one would run into a problem. 
In this limit, the vortex phase drops out as a dynamical 
degree of freedom, and the photons are counted in the 
same way as in the Coulomb phase except that they are 
now massive: the interactions between static vortices are 
now screened, while the transversal 'spin wave' photon 
acquires a mass. In this way one would find one propa- 
gating massive mode instead of two. 

In the relativistic theory (cy — c) it fits neatly: we 
encounter now four fields {At , A^ , Ay , (py) subjected to 
one gauge constraint. Of the remaining three physical 
fields one takes care of static interactions, and we are 
left with two photons which are degenerate! Since only 
one of these photons is carried by the smooth phase field 
configurations, the other one reflects the phase rigidity 
of the vortex condensate: one of the two excitons light- 
ing up in the phase world is just telling that the dual 
superconductor resists a twist in its phase. 

Having given away the bottomline, let us now substan- 
tiate these matters with some explicit computations. 

III. DISORDER FIELDS PROBED BY ORDER 
MEANS 

It does not seem to be widely recognized that the prop- 
agators of order fields can be straightforwardly expressed 
in terms of the disorder fields. This was discussed in 
an earlier paper dealing with field-theoretical elasticity^, 
and to save the reader the effort of learning this intricate 
affair, let us rederive these relations for the far simpler 
Abclian-Higgs case. 

These relations are associated with the first step 
of the duality, where the phase fields are turned 
into photons mediating interactions between the vor- 
tices. To set the stage, let us shortly review these 
mattera^^l^^^l^°1"l^^1"l^^1^^l^^. After coarse graining, 
phase dynamics can be written in terms of the La- 
grangian, 

The coupling constant g is proportional to the original 
coupling constant g. The spin-wave velocity is given by 
ratio of the stiffness and compression moduli (? = Ps/ks 
and set to 1 in the last step. It is left implicit that the 
phase field is compact, = 0-|-27r. We take the superfluid 
'three' velocity 

^;^(x) = 9^0(x) (4) 

as the natural, 'primitive' observable of the orderly side. 
In the ordered phase, the momentum space velocity- 
velocity propagator is proportional to the phase-phase 
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propagator, {{Vf^\v^))q,^ = ql5^^ {{(t)\'t>))q,Lj and the lat- 
ter sufRces to calculate the order parameter propagator 
((e"^|e"^)) (e.g., ref.'^). In the disordered phase (j) itself 
becomes multi- valued and meaningless, but continues 
to be single valued and meaningful. 

In the phase-ordered state the theory Eq. is Gaus- 
sian and the velocity propagator is easily computed by 
adding a generating functional to the Lagrangian, 



followed by taking the functional derivative 
1 d^Z[J^] 



(5) 



(6) 



The non-relativistic propagator measured in condensed 
matter experiments represents only the subset of compo- 
nents of the relativistic propagator Eq. (jSJ with spatial 
indices: {{vi\vj)). In the phase ordered state of the XY 
model one can integrate the Gaussian Goldstone fields in 
Eq. (O with the result. 



z[j,] = n 



2TTg 



- 7 - 



(7) 



and the propagators follows immediately from Eq. Q 
identity. The relativistic and non-relativistic versions are 
respectively, 



{{vt,K)) = g 



2 I „2„2 «J ■ 



(8) 
(9) 



where the longitudinal- and transversal (for later use) 
projection operators are. 



pi 



«2 ■ 



„2 ' 



(10) 



Obviously, we find only a single mode: the single Gold- 
stone boson of the scalar theory as shown response plot 
Fig. 

Of course, the above procedure no longer works in the 
absence of the phase condensate. At any finite disorder 
there are configurations present containing topological 
defects (vortices) , which are ignored in the path integral 
Eq. lO but these are easily handled in the language of 
the dual disorder field theory. 



Let us turn to the duality itself. The first crucial step is 
a simple Legendre transformation. Introduce Hubbard- 
Stratanovich auxiliary fields such that the phase ac- 
tion including the external sources, Eq.(|3J) is written as. 



f Jm^m + iaM^. + f + ^^^^^^4> (11) 



Let us ignore the external sources J^^ for a moment, to 
focus on the duality itself. The phase field is split into 
smooth- and multivalued pieces: = (jigm + 4>mv- By 
reshuffling derivatives {^fid^(t)sm = -4>smd^ipi,) the Gaus- 
sian (psm turn into Lagrange multipliers imposing a con- 
servation law on the auxiliary fields. 



(12) 



In the phase dynamics interpretation, just represents 
the supercurrents and Eq. (|12|l is the hydrodynamical 
continuity equation governing superflow. In 2+ ID conti- 
nuity can be imposed on the supercurrents by writing 
in terms of non-compact C/(l) gauge fields as 



(13) 



and it follows that the remaining pieces of the Lagrangian 
Eq. (|ll|l (ignoring the J'^s) can be written as, 

+ ^^Mmv = ^^F^.F^- + lA.Jl (14) 

where F^^i, = d^A^, — di,A^ is the field strength of the 
dual gauge sector while 



J.. 



d^dx(f)MV 



(15) 



is the non-integrability of the phase field, having the 
meaning of vortex current. This completes the key step 
of the duality: it shows that in 2-l-lD the rigidity of the 
phase medium can be exactly parameterized in 'force car- 
rying photons' (the A^'s) while the vortices act as sources 
for these photons. The disorder field theory Eq. Q 
follows immediately: vortices are indistinguishable from 
charged bosons interacting with electromagnetic fields 
and a system of such bosons is described by the Ginzburg- 
Landau- Wilson theory"*. 

This is all familiar territory but the following simple 
operation seems not to be commonly known. Let us in- 
clude the external sources and use the identity Eq. to 
calculate the phase propagator, but now using the action 
after the Hubbard-Stratanovich transformation, 



Ju=0 



9^ f dx^ {^J^J^-igJy.i^-^i^iy.-ii^dy. 



(16) 
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This implies an exact relationship between the veloc- 
ity propagator and the propagator of the supercurrents, 
rooted in the Legendre transformation, 



(17) 



Because = e^j^^xduAx this implies that in fact the 
phase velocity/spin wave propagator is proportional to 
a linear combination of the physical photon propagators 
of the dual gauge disorder-field theory. This implies that 
the poles of the magnon and photon propagator have 
to coincide and this has to be because both describe the 
same physics. However, the pole strengths might be quite 
different reflecting the 'dual relativity principle': pend- 
ing the use of order or disorder 'tools' one might get a 
very different view of the same underlying reality. The 
result Ea. (|17|l . first derived in Ref.'', shows that at least 
in the Abelian-Higgs case, the two observers should ac- 
tually agree more on what they see than one could have 
expected a-priori. The key is that although the vortex 
condensate falls prey to dual censorship, the orderly ob- 
server can still learn much about the dual world because 
he/she can probe the dual photons according to Ea. (fT7|l . 

IV. MAGNONS AS PHOTONS 

Let us exercise the notions of the previous section in 
the simple case of the phase ordered state. We know 
the answer (the Goldstone mode, Eq. (S)), we know the 
dual side (Maxwell theory, Ff^^.F^'^), and we know how 
these relate (the 'Zaanen-Mukhin' relation, Eq. I17|) . ft 
is indeed a straightforward exercise. 

Although the dynamics is fully relativistic the ques- 
tions of relevance to condensed matter experimentalists 
are not relativistic: only the spatial components of the 
propagator are measurable (Eq. ^ . This makes it conve- 
nient to use a Coulomb gauge fix. The Maxwell action in 
momentum-Matsubara frequency space becomes, includ- 
ing the external sources J^^*', 



EM 



1(4, At) 



(j2 -icj„(q| 
iJ''^*- • At. 



A 



(18) 

where we have explicitly indicated the time- (A^t) and 
space (X) components of the gauge fields and currents. 
The bra- and ket in the gauge field propagator represent 
rows and columns g.i, respectively. Notice that from now 
on we keep the Goldstone/spin-wave velocity c explicit, 
for purposes which will become clear later. 

Provided that we choose a gauge fix T that does not act 
on the temporal component Ar , the temporal component 
can be integrated out. This yields the usual Lagrangian 
with Coulomb interactions between static sources. 



^EM — 



/t 7 
+i{JL - 



Atp^A 



(19) 



The longitudinal component A^ is unphysical (its source 
is iujnJr — qJh — > OtJt + diJi = 0) and it should be 
removed by the Coulomb gauge 



= d,Ai = qAL. 



(20) 



We end up with two propagators for the gauge fields, as 
it should in 2-l-lD. We find one dynamical photon. 



((4 14)) 



(21) 



and a propagator taking care of the Coulomb interactions 
between the static sources. 



((4l^r)) 



1 

gq^ 



(22) 



This is of course textbook electromagnetism, but be 
aware of the twist in the interpretation. The photons 
now keep track of the capacity of the phase condensate 
to respond to external influences. The outcome is: it is 
carrying a 'Goldstone photon' {At, Eq. 121(1 and it can 
mediate as well interactions between static vortices {At, 
Eq. 1^ . 

We are now in the position to evaluate the 'Zaanen- 
Mukhin' relation Eq. H17|l . For this purpose, we are only 
interested in the spatial components of the supercurrents 
The supercurrent propagator is easily found by using 
the definition Eq. H13|l . and the results for the gauge 
field propagators Eqs. 1)21112211 and wc find for its spatial 
components. 



lot. 



2 y 



Pi 



(23) 



Using now the Zaanen-Mukhin relation Eq. (|17|l . 



n pL 

„2 2 ij 



P^ 



g 



9 9 

c q 



uj^ + c/q 



,2„2 U 



(24) 



After this long detour, we indeed have managed to re- 
cover the spin wave propagator Eq. ^\ 

The simple lesson following from this simple exercise 
is that the dual photon language is in a way more com- 
plete than the description in terms of phase fields, in the 
sense that the gauge fields keep track in an explicit way 
of both the capacity of the medium to propagate Gold- 
stone bosons and the fact that it mediates interactions 
between its topological excitations. The Zaanen-Mukhin 
relation filters out the Goldstone sector from the 'om- 
nipotent' dual gauge sector, keeping its topological side 
(the Coulomb propagator, requiring vortex sources) com- 
pletely hidden from the eye from the 'orderly' observer. 
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V. WATCHING THE DISORDERED STATE 
WITH ORDERLY MEANS 

Surely, the dual route of the previous section is a rather 
inefhcient way to derive the propagator of a Goldstone 
mode. This changes drastically in the phase disordered 
state. Resting on the fact that the dual gauge theory is 
now governed by order, precise information on the second 
sound propagator can be extracted with barely any extra 
investments. The only other option is the strong coupling 
expansion in the Hamiltonian language and this becomes 
very tedious at intermediate couplings. 

The phase disordered state corresponds with the Higgs 
phase of the gauge theory Eq. Q , corresponding with the 
state where vortex loops have blown out and the vortices 
have Bose-condensed. As a consequence, the bosonic dis- 
order field ^' = l^ole**^^ acquires a finite expectation 
value. This theory is fully relativistic, as we explained, 
and this vortex condensate is literally like the U{1) Higgs 
phase of high energy physics^^. It will turn out to be 
quite convenient for the interpretation of the results to 
consider a non-rclativistic extension of the theory char- 
acterized by a condensate velocity cy, which is different 



from the spin-wave velocity c, entering the time compo- 
nents of the covariant derivatives ~ -^{9^- — iAt). Of 
course, the real theory is characterized by cy = c as im- 
plied by the Lorentz-invariance of the action Eq. (|3J). 

Let us employ the usual unitary gauge, corresponding 
with fixing the condensate phase 4>v = 0- The finite ex- 
pectation value of the disorder field results in the familiar 
Higgs term in the action 



■^1 



Higgs 



(25) 



The only specialty is the velocity cy. In high en- 
ergy physics this is the light velocity while in the non- 
relativistic condensates of condensed matter physics cy 
is the sound velocity (in BCS theory ~ iipiSiiS), which 
is vanishingly small compared to the light velocity with 
the consequence that one can get away with a time inde- 
pendent Ginzburg-Landau theory. In our final result we 
have to set cy = c. 

We now follow the same route as in the previous sec- 
tion. Adding the Higgs term the Lagrangian becomes. 



- f(4,At)^' + f .J^;^2K■2pT)(^)+^^rAl + ^J■Al (26) 



Introducing a Higgs mass f2, defined by 



(27) 



Since the gauge has already been fixed, the temporal 
components At can be safely integrated out. 



= f At 



r!2) 



'^yq 



r!2 



At 



JlJr 



2g 2 



.(28) 



The last term corresponds with the interactions be- 
tween the static vortices which are now short ranged. 
The interest is in the dynamics of the gauge fields itself. 
As before, we find a transversal photon At characterized 
by second sound propagator which has acquired a Higgs 
mass. In addition, we find an extra longitudinal photon 
(the first term) which is now physical. This is also char- 
acterized by the same Higgs mass but it is propagating 
at the condensate velocity showing that it represents the 
phase rigidity of the dual superconducting matter sector. 

The propagators for gauge fields are easily determined 
from the inverse of the full action Eq. (|26|) . The super- 
fiuid current propagator is decomposed into longitudinal 



and transversal parts ^^'^ (parallel and perpendicular to 
the momentum q respectively) and the propagators are 
found to be 



((aia)) 

{{£,T\iT)) 



c^q^ 



(29) 
(30) 



Using now the Zaanen-Mukhin relation Eq. H17|) and 
the momentum propagators Eqs. H29I - 1^0 ) 1 . we obtain 
the result for the non-relativistic propagator for the su- 
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FIG. 2: Spectral functions associated with the superfluid 
velocity-propagator, computed on the Gaussian level using 
the dual theory . These results should be accurate deep in- 
side the ordered and disordered phase (a) The ordered (su- 
perfluid) phase: the second sound pole of Fig. ^ is recovered 
as it should, (b) The disordered (Mott insulating) phase: we 
used here a condensate velocity (cv) which is half the sound 
velocity c for the mere purpose to make visible the different 
behaviors of the strength of the second sound- (higher branch) 
and dual (vortex) condensate (lower branch) poles. In reality 
these velocities are the same and the modes are degenerate. 
For q — + the pole strengths of the two modes are the same, 
while they are governed by the same Higgs mass, and they 
can be combined in the ±1 helicity modes as expected from 
the strong coupling expansion in the Hamiltonian formalism 
(Fig. Id; see section V). However, for increasing momen- 
tum the condensate pole looses gradually strength while the 
second sound pole becomes more and more like the 'orderly' 
result of (a), reflecting that at distances short compared to 
the dual London penetration depth the medium 'rediscovers' 
the order. 



perfluid velocity in the disordered phase. 



(31) 



The spectral response from this propagator is plotted on 
Fig. Eb. 

The longitudinal (first) term represents, as before ( Eq. 
I24|l . the correlations associated with the smooth part of 
the phase field : this is literally second sound acquir- 
ing a mass associated with the disappearance of the su- 
perfluid rigidity at large lengths- and times. We notice 
that in the static limit (w„ — > 0) the longitudinal part 
becomes a constant, signaling that even at the shortest 
distances superfluid correlations have disappeared. This 
makes sense: when vortices populate the whole system, 
then any long living correlation is destroyed even between 
two neighboring sites when one waits long enough. 

The second, transversal term is the interesting one: we 
indeed find a second mode and although it has the same 
mass as the gapped second sound it propagates with the 
condensate velocity. It is of course the longitudinal pho- 
ton reflecting the dynamics of the dual superconducting 
vortex matter. In order for the superfluid velocity cor- 
relator to acquire a non-zero transversal component it 



is actually a requirement that the phase fleld becomes 
non-integrable. This becomes clear by inspecting the 
transversal part of the supercurrent Eq. (^0)1 , 
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m + di(f>Mv) = - 



(32) 



>MV 



where the smooth part has disappeared since it makes no 
sense to have derivatives in the transversal direction for 
the case of smooth flelds. 

Let us analyze the result Eq. H31|) in more detail. 
The velocity cy has done its job in establishing that 
the transversal poles of Eq. (|31|l are indeed due to 
the superconducting vortex matter, and we can impose 
Lorentz invariancc by setting cy = c: the longitudinal- 
and transversal modes become degenerate as they should 
because of the degeneracy associated with the n = ±1 ex- 
citons, trivially seen in the Hamiltonian formalism. A dif- 
ferent issue is the pole strength of the second sound (J^) 
and condensate poles {It) as measured by the velocity- 
velocity correlator, characterized by the ratio. 



II 
It 



(33) 



Giving this a minute of thought this makes a perfect 
sense. First, aX q — Q both modes at real frequency 
Lo = 0, have the same weight. It follows from the Hamil- 
tonian formalism that at wavelengths large compared to 
the London length the excitons correspond with the he- 
licity ±1 eigenstates of the angular momentum of the 
0(2) quantum rotors. The supercurrents have the status 
of canonical momenta and should therefore be combined 
in currents with definite helicity ~ ± iS^t- The 
implications is obvious: at q — > the longitudinal- and 
transversal poles of the velocity propagator should have 
the same strength because all what exist in this limit are 
the helicity eigenstates. 

What is changing at smaller distances? The character- 
istic momentum scale is of course the inverse dual London 
penetration depth q^ ^ — ^/c and from Eq. 

it follows that at larger momenta the strength of the dual 
condensate pole decreases quadratically in momentum 
relative to that of the second sound pole. Within the con- 
fines of this gaussian treatment this makes again sense. 
At these short times and distances one enters a regime 
where one is probing mainly the phase ordered matter 
forming the background in which the vortices move. This 
matter is the same stuff as the fully phase ordered mat- 
ter and accordingly it should carry the same Goldstone 
excitation. Eq. tells how to interpolate between the 
disorder physics at (7 — )■ with its number eigenstates 
and the phase ordered regime at large momenta: the 
dual condensate pole looses its weight gradually, in fact 
in the same way as the Higgs mass looses its influence on 
the dispersion. 
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VI. DUAL VIEW ON THE CRITICAL REGIME 

We are not done yet. We have implicitly assumed up 
to this point that the fields are non-interacting. Mod- 
ulo perturbative corrections this would have been fine in 
dimensions above the upper critical dimension but the 
Abelian-Higgs model in 2+lD is below its upper critical 
dimension due — 3+ ID. One has now to be cautious with 
considerations like the one in the previous section. Upon 
Exceeding the scale il one does not re-enter the ordered 
phase but instead one enters the quantum critical regime 
which has no longer to do with order or disorder but has 
acquired an own identity due to the strongly interact- 
ing nature of the critical point. Away from the critical 
coupling, the order- (second sound) and dual order (lon- 
gitudinal and transversal photons) excitations discussed 
so far still make sense because they will appear as bound 
states pulled out from the low energy side of the contin- 
uum of critical modes, with a pole-strength and binding 
energy diminishing upon approaching the critical cou- 
pling. The missing link at this point is the appearance of 
the continua of critical modes as picked up by the veloc- 
ity propagator. After some preliminaries we will derive 
their form resting on the large body of knowledge on the 
3D XY critical state. These critical continua turn out to 
behave in a quite surprising way, with the second sound- 
and condensate contributions showing a completely dif- 
ferent behavior away from q = (see Fig. Oand Fig. 
We will subsequently focus in on the detailed way the 
quasiparticle poles (second sound, the excitons) develop 
as function of the distance from the critical point, mak- 
ing the case that the critical continua have to be as they 
appear in order to be consistent with the quasiparticle 
poles. 

In order to describe the system close to- and right 
at criticality, we introduce renormalized parameters and 
critical exponents. The role of reduced temperature is 
taken by the quantity e = which is the reduced 

coupling constant. If e < or e > 0, we approach the 
quantum critical point at gc from the ordered and dis- 
ordered side, respectively. Pending if we approach the 
critical point from the order- or disorder side, the system 
will scale either to the stable fixed points associated with 
phase order and non-interacting second sound {g = 0) 
or with non-interacting rotors {g = oo). The reduced 
coupling constant e is therefore a relevant operator with 
scaling dimension i/f^ > 0. Another relevant field, which 
plays the role of the magnetic field in the standard scal- 
ing analysis is the generating functional field . Since it 
is relevant at the transition, its scaling dimension is also 
positive yj > 0. 

The model we consider is relativistic, with dynamical 
critical exponent z = 1, and its critical behaviour will co- 
incide with that of the 3D XY-model. We use the state 
of the art for the exponents, based on analytic methods 
(high-temperature expansionSS, vortex-loop scaling^i, 
one-loop renormalization groupSS) as well as numeri- 
cal results from Monte Carlo simulationsSiSi24. The 



critical exponent rj for the order parameter propagator 
(e^'^^e-*'*') has been studied in great deta.i^i^i2^i2ML2^. 
However, our interest is in the velocity correlation func- 
tion Eq. which is not straightforwardly related to the 
vertex correlator away. Instead, we will use the knowl- 
edge of the scaling dimensions of the dual field r]A to 
derive the form of the velocity propagator in the critical 
regime. 

Let us first analyse the model and its propagators right 
at the critical point g = gc- The exponent rjA is usu- 
ally defined as the critical exponent of the gauge fields 
correlation function right at the critical point g = gc, 
i.e. {{AA}) (X l/p^"''-*. To be consistent with the 
literatur o^i^^ , we have to change the gauge fix from 
'our' Coulomb/unitary gauge fix to the Lorentz gauge fix 
(d^Af^ — 0, i.e. the vector potential is purely transver- 
sal). In this gauge fix, the gauge field can be projected 
onto a 3D linearly polarized basis (defined as Gq — ^, 
e_i = —bt and e^-i = e_i x Bq). The component 
is set to zero by the gauge fix, with only the space- 
time transversal components of the fields being physical. 
The spatially transversal photon (second sound) degree 
of freedom At is now represented by A^i. The remain- 
ing component A+i that admixes the Coulomb and the 
longitudinal photons plays the role of the vortex phase 
degree of freedom in this particular gauge fix. On the 
Gaussian level of the previous sections, the propagators 
for the gauge fields within the Lorentz gauge fix degen- 
erate and are given by 

,Jl'''o2 - (34) 

In the Coulomb phase one finds the same propagator with 
a vanishing gap, il = 0. The indices are taking 'transver- 
sal' values h,h' = ±1. The coupling constant in the 
prefactor is expressed in terms of the superfluid stiffness 
Ps = 1/.9 which is a quantity which does renormalize. 
It follows that the residues of the quasi-particle poles 
(order/disorder excitations) are also renormalized which 
would not be the case if the prefactor would correspond 
with the bare coupling g. The overall prefactor in the 
expression for the velocity propagators corresponds with 
g^Ps in this scheme. The g^ is the bare critical coupling 
since the relation between the dual and original propa- 
gators Eq. I|17|l is an exact relation from the Legendre 
transformation, which is also valid in the critical regime. 
Accordingly, both the second sound of the ordered side 
and the excitons of the disordered side loose their pole 
strength approaching the critical point and this is gov- 
erned by the renormalization of the superfluid density ps 
which we will deduce starting from the known critical 
behaviour of the dual gauge field propagators. 

Herbut and Tesanovio^^ analyzed the charged XY 
model which is equivalent to the dual action Eq. jSJ. 
From their expression for the /3-function governing the 
renormalization of the electrical charge, it follows that at 
the fixed point 

= eliD-A + Ju)- (35) 
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Assuming that the charge scales to a finite value eg, it 
follows that rjji = 4 — D = 1. The same result was 
obtained by Hove and Sudb0^, using Monte-Carlo to de- 
termine the exponent r]A from the vortex correlations at 
the critical point. They introduced a relation between 
the correlation function of the vortex tangle G{p) and 
the dual gauge field propagator, 



(AtA) 



2/3 2(37r^Gip) 

,2 \ ~9 



p. 



(36) 



valid for the case of the uncharged original/charged 
dual action. Notice that in Refi^ h is used for the 
dual gauge fields and A for the original gauge fields. 
At the critical point the vortex correlator is given by 
limp^o 2(3-k'^G{p) = p^ - C3{g)p'^+'^^ + ... and using Eq. 
(EEll it follows that 



p'{A^A)^C3{g)p^ 



(37) 



According to their numerical simulation, Eq. 1)3 7|l shows 
a linear behaviour and Hove and Sudb0^ conclude that 
VA = 1- 

The critical propagator of the dual gauge fields, Eq. 
(|37|l . can be used to establish the form of the velocity 
propagator Eq. ^ in the critical regime. Comparing 
Eq. H37|l with our form of the gauge field propagator Eq. 
(|34|l . and bearing in mind the degeneracy, we conclude 
that each field component propagator corresponds with 
one half of the propagator Eq. 1)3 7|l . 



-P 



(38) 



I 



7? 



ql<q2<q3<q4 



CO 



FIG. 3: The spectral functions associated with the second 
sound (full lines) and condensate (dashed lines) pieces of the 
velocity propagator in the critical regime for various momenta 
(Eq. MU l. At g = both critical continua become degenerate 
and linear in to reflecting the simple correlation function expo- 
nent riA = 1 associated with the dual gauge fields. However, 
at finite momenta it is seen that the second sound continuum 
diverges at threshold while the condensate piece is actually 
suppressed, although both continue to be governed by the 
same scaling dimension. At finite momenta this different be- 
havior of the critical continua has to be present in order for 
them to be consistent with the momentum dependence of the 
pole strengths of the propagating disorder excitations appear- 
ing at the moment one moves from the critical coupling. 



We can now use again the universal Zaanen-Mukhin re- 
lation Eq. (|17|l to obtain the velocity propagator 



2-riA 



P- 



(39) 



right at g = gc- This is the first main result of this sec- 
tion. The dots represent constant terms with no imagi- 
nary parts as well as short distance corrections. At least 
deep in the critical regime the Wick rotation to real time 
is simple- because scale invariance implies that Euclidean 
propagators are power laws, turning into branch cuts in 
real frequency With 77^1 = 1, right at the criticality, the 
spectral function has two quite different branch cuts in 
the longitudinal and the transversal channel 



lm{{vi\v.j))T ^ 0{u. 



„2„2\ 



,2„2 



(40) 
(41) 



and we sketch both pieces of the velocity correlator in 
Fig. (jSJ. 0{x) is the Heaviside unit step function. 

This is quite an unexpected result. At q' = we find 
both spectral functions to be simply proportional to fre- 
quency, a simple behavior which of course originates in 



rjA — ^- Upon increasing momentum, the 'sound' and 
'condensate' spectral functions start to behave very dif- 
ferently near the threshold uj = cq although at large uj 
they merge together again. The sound part shows the 
usual& divergence uP' jiio — cq)'^~'^^ ~ uP jiuj — cq) while 
the condensate piece develops like (w^ — Pq^)'^^/'^ = 
^ up — . The degeneracy of the two contributions at 
q — Q rings a bell: at infinite wavelength it should be that 
the critical fluctuations are eigenstates of rotor angular 
momentum, 'equalizing' the condensate and sound con- 
tributions as we found for the propagating excitations. 
To understand better why these contributions should be- 
come different at finite momenta, we should first analyze 
in more detail what happens with the quasiparticle poles 
close to the critical point. 

To analyze the behavior of the quasiparticle poles in 
the ordered- and disordered phase close to the critical 
coupling, we need hyperscaling. Although one has to be 
careful^—, recent numerical simulations^^ show that there 
is none or a very small violation of the hyperscaling for 
3D XY. Let us first repeat the standard hyperscaling ar- 
guments applied to the velocity-velocity propagators. We 
denote the propagators of the gauge field in real space as 
6*^1(0;, e). It is generated by the term JhJ^h in the action 
and this generating functional of the gauge fields Jh plays 
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a role similar to a magnetic field. It is a relevant field 
that scaling dimension yj. Hyperscaling requires that 
such fields act on a block of 6'^"'"^ points in space-time, 
treated as a single variable. After a scale transformation, 
the new propagator is related to the original one by 



L2(d+1) 

—^G{r,J) (42) 



A 



J 



where Xj = e^^ is the scaling factor of the generating 
functional for the gauge fields. 

Repeating the scale transformation n times in the 
vicinity of the critical point, we obtain 



\2n 



(43) 



with the scaling factor — e^' associated with the re- 
duced coupling constant. Choosing n such that (Ae)"5 — 
const, from Eq. (|43l) it follows that the propagator be- 
haves universally on both sides of the critical point as, 

GA(r,e) cx \e\^'-'^+^-y^'^^±{r/\e\-^). (44) 

Te Functions $± are universal functions associated with 
the ordered and disordered sides of the critical regime, 
and given in terms of G{r, const). The denominator in 
its argument is the correlation length that diverges at 
the critical point with exponent v, implying the familiar 
relation ly = I /ye- The relation of the 'magnetic field' 
exponent yj to the scaling exponent rjA follows when we 
set e = in Eq. (O, 



yj 



d + 3 - T]A 



(45) 



using the known value rjA = I. Together with the relation 
for ly, Eq. (|44() can be written as 

GA{r,e) oc |e|'^(''-i+''-*)$±(r/|e|-'^) (46) 
- kP''$±(r/|er^). (47) 

This is just the familiar result that the behaviour of the 
correlation function close to the critical point is governed 
by the exponents v and 77 (with 77 = 1 in the present case), 
and the crossover functions $±. 

Let us first approach the critical point from the dis- 
ordered side, i.e. e — > 0"*". This phase is characterized 
by the gap Eq. H27|l . which we can call (compare Ref.^) 
A_|_ = Q. This gap is proportional to the inverse correla- 
tion length of the vortex tangle ^ = c/fl. Upon approach- 
ing the critical point, both the correlation length vanishes 
and the gap diverges with characteristic exponent v as 
^ oc and A_|_ cx e^", where z = 1 is the dynamical ex- 
ponent which equals one in this specific case. The 3D XY 
correlation length exponent 1/ = 0.66 — 0.67 « | accord- 
ing to a large body of workSfliSiiSiSi. Given that there 
are two dynamical fields in the problem {A^ and the vor- 
tex phase field phase 4>v) one could be tempted to think 
that there are two correlation lengths in the problem, 



but this is not the case the problem is effectively Lorentz 
invariant, consistent with the numerical worlsiSiSS*^^ and 
an argumenliSi linking it to the anomalous dimension of 
the gauge field rjj^i^. 

The scaling dimension of the superfluid density can be 
deduced from Eq. (|34l) . After Fourier transformation to 
space-time, the Gaussian propagator Eq. behaves 
like 



(48) 

p,r<'^'^*+(|) (49) 
p.e''('^'i)$+(|). (50) 



Comparing it with the hyperscaling form for the gauge 
field propagator Eq. 144|l we conclude that the superfluid 
density scales as 



Ps oc 



|i>(2-77a) 



(51) 



at the disordered side of the critical point. 

We have now arrived at a point where we can de- 
termine the behaviour of the two quasi-particle poles 
upon approaching the critical point from the disordered 
side. Using Eq. H39|l . the fact that g g^ps and the 
scaling of both ps and we conclude that the vortex- 
condensate pole ^ has a strength proportional to 
Ps^^ oc 62ziy+iy(2-j)A) _^ g3i^^ vanishing upon approach- 
ing the critical point with an exponent = 2 while its 
strength disappears in the critical continuum as indicated 
in Fig. Turning now to the second sound pole P^, 
we observe that at long wavelength {q 0) its strength 
behaves exactly like the condensate pole. This has to be 
because eventually, at large enough distances, one should 
recover the fact that these excitons correspond with the 
exact rotor angular momentum eigenstates. However, for 
increasing momenta the term in the numerator ~ c^q^ 
takes over, and the strength of the large momentum sec- 
ond sound pole is scaling more slowly to zero upon ap- 
proaching the critical point, governed now by the super- 
fluid density exponent ^{2 — tja) = |. This is of course 
not different from what we found on the Gaussian level, 
with the second sound pole overtaking the condensate 
pole when the vortex condensate is 'loosing its grib', gov- 
erned by the Higgs mass ft. 

To complete the picture, let us finally consider what 
happens with the second sound pole approaching the crit- 
ical point from the ordered side. This is straightforward: 
as before, we should substitute g ^ g^Ps in the gaussian 
result Eq. Q and ps le]" because ps renormalizes in 
the same way on both sides of the transition^S. In other 
words, the strength of the second sound pole on the or- 
dered side coincides with its behavior at large momenta 
on the disordered side. 

Not surprisingly, we have found that the 'order poles' 
behave quite like the results we found on the Gaussian 
level in the previous sections except that renormalized 
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mass scales and -quasiparticle residues have to be used, 
all governed by the same correlation exponent v , because 
riA 'magically' drops out. We can now use this knowledge 
to comprehend why the critical continua of Fig. |31behave 
the way they do. We already argued that at energies far 
away from the threshold lo = cq the second sound and 
vortex condensate pieces picked up by the velocity cor- 
relator merge in the same linear I lo behavior. At 
finite q the differences between the two are large near 
threshold. With the knowledge regarding the behavior 
of the quasiparticle poles at hand this now makes sense. 
Ps being a relevant operator, its influence at high ener- 
gies is small while growing when times get longer. A bit 
away from the critical point, it takes over at a length 
~ ^ where the system gets under control of the stable 
fixed points at zero or infinite coupling, which are also 
in charge of protecting the quasiparticle poles. Surely, 
the quasiparticles close the critical point can be viewed 
as bound states pulled out of the critical continuum due 
to the effect of the relevant operators (Fig. How- 
ever, because of the way the latter scale, the quasipar- 
ticles are formed from the low energy end of the critical 
continuum. What does this mean for our velocity prop- 
agator, 'watching' the true critical excitations through 
the 'duality filter' ? We derived some clear rules for how 
the weights should be distributed over the quasiparticles: 
the condensate and sound poles of the disordered state 
should have equal weight at q = 0, but the former should 
loose its weight rapidly for increasing momentum. In- 
specting now the low energy end of the critical continua 
for various momenta we see this rule also at work (Fig. 
O! We notice that this 'weight-matching' of the criti- 
cal continua and the quasiparticle poles is to an extent 
even quantitative. For this purpose we inspect the pole 
strength ratio Eq. close to the critical point. For 

fixed q, due to the gap in the denominator, we learned 
already that the ratio diverges like ^ q^ je^^ . However, 
the prefactor of the second sound pole strength is how- 
ever proportional to . Comparing it with the ratio of 
the spectral responses right at the critical point and near 
threshold (w « cq) 

(t^ = l' X "divergent part". (52) 

\It J uj^~c^q^ 

We find a perfect match - the strengths of the spin- wave 
and the condensate excitations are proportional to the 
strengths of their respective critical continua where they 
have their 'origin'. 

Surely, this does not explain everything, and to a de- 
gree Eq. (|4HI is a result which stands on its own. How- 
ever, given the simple integer 7]a exponent, it appears to 
us to be a unique analytical form which obeys the general 
requirements of scale invariance and Wick rotation, hav- 
ing at the same time the right form to be consistent with 
the evolution of the spectral weights in the quasiparticle 
poles. 
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FIG. 4: Cartoon of the appearance of the second sound (full 
lines) and dual condensate (dashed lines) contributions to the 
velocity spectral functions at finite momentum in the close 
vicinity of the quantum critical point, both on the ordered- 
(left) and disordered (right) side. Although the critical con- 
tinua are expected to have a very similar appearance, on the 
ordered side only a second sound pole is found. However, on 
the disordered side the system scales to dual superconducting 
order with the effect that one finds both a propagating second 
sound- and condensate excitations with strengths governed by 
the XY correlation length exponent. However, the way their 
pole strengths develops as function of momentum tracks the 
gaussian result shown in Fig. 2b and it turns out that the 
momentum dependence of the critical continua is just of the 
right kind to be consistent with the behavior of the disorder 
poles (section VI). 



VII. CONCLUSION 

What have we learned? We have taken the simple ex- 
ample of phase dynamics at zero chemical potential to 
predict the form of the superfluid velocity correlator in 
the ordered, disordered and critical regimes, exploiting 
the vortex duality. Although it does not seem to be gen- 
erally recognized, the dual order of the disordered side 
does manifest itself when the system is interrogated with 
'orderly means'. In our example, one of the two degen- 
erate excitons of the phase disordered/Mott insulating 
state can as well be called the longitudinal photon asso- 
ciated with the phase-rigidity of the dual superconduc- 
tor. Although less obvious in the present simple exam- 
ple, this might be used as a technical convenience. Quite 
generally, it is easier to compute the excitation spectrum 
of the system 'around' the ordered state, helped by the 
Goldstone theorem, the Higgs mechanism, etcetera, while 
disorder is not a convenient starting point when the in- 
teractions are strong, ft is actually so that this work was 
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originally inspired on problems encountered in the study 
of quantum liquid crystals where this 'dual route' to the 
spectral functions associated with measurable quantities 
seems to be the only way available''. 

This dual route also tells another story which is far less 
obvious. This can be summarized as: 'studying the disor- 
dered state with order operators, one recovers the signal 
of the ordered state at energies and momenta where the 
dual order parameter looses control.' In our specific ex- 
ample, the dual condensate piece of the exciton doublet 
of the disordered state fades away when momentum is 
increased and at large momentum only second sound re- 
mains. We argued that this same 'mechanism' is even 
at work in the critical state, being ultimately responsi- 
ble for the rather odd appearance of the continuum of 
critical modes as measured by the velocity correlator. 

Within the field-theory this is surely correct - is is 
based on controlled calculations. Another issue is, can 
we literally apply the field theory to condensed matter 
problems for this particular purpose? We are actually 
not sure. One way to read the effect of the previous 
paragraph is as follows. It is assumed in the field theory 
that the lattice constant is vanishing. As applied to any 
problem with a finite lattice cut-off, this means that the 
field theory can only be taken literally when the distance 
between vortices is large compared to the lattice constant 



- the small fugacity limit is implicitly taken. If this is 
the case, by zooming in one will eventually get at length 
scales which are smaller than the inter-vortex distance, 
and here one will rediscover the implacable order and its 
dynamical implication in the form of its Goldstone mode. 
The way this limit is reached is a bit more sophisticated 
than suggested with these words, but this we discussed 
at length in this paper. 

The problem is that in real condensed matter systems 
the vortex (or, in general, 'dual') matter is actually rather 
dense when density is measured in units of the lattice 
constant. Although the field theory is fine in the scaling 
limit, as helped by universality, to what extent can the 
effects we discussed in this paper become noticeable? We 
leave this as an open question, in the hope that others 
might get motivated to have a closer look. 
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